Observation of a subgap density of states in superconductor-normal metal bilayers 

in the Cooper limit 
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We present transport and tunneling measurements of Pb-Ag bilayers with thicknesses, dpb and 
dAg, that are much less than the superconducting coherence length. The transition temperature, 
Tc, and energy gap. A, in the tunneling Density of States (DOS) decrease exponentially with dAg 
at fixed dpb. Simultaneously, a DOS that increases linearly from the Fermi energy grows and 
fills nearly 40% of the gap for Tc ~Q.l T^^uik- This behavior suggests that a growing fraction of 
quasiparticles decouple from the superconductor as — >0. The linear dependence is consistent 
with the quasiparticles becoming trapped on integrable trajectories in the metal layer. 
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Simple metallic phases have a finite resistance in the 
zero temperature limit and a nonzero Density Of elec- 
tronic States (DOS) at the Fermi energy, Ep. In two 
dimensions (2D), however, the scaling theory of local- 
ization asserts that simple metallic phases do not exist 
. Ultrathin films exhibit a continuously decreasing con- 
ductance with decreasing temperature in support of this 
assertion Q . Nevertheless, there exist an increasing num- 
ber of quasi-2D systems with metallic transport proper- 
ties at low temperatures. These include 2D electron gases 
formed in semiconductor heterostructures, which appear 
to show an Insulator to Metal transition with electron 
density ^] and ultrathin films of metals balanced on the 
brink of a superconducting transition by disorder j^j or 
magnetic field These metallic phases are probably 

not simple Fermi liquids and their existence depends on 
electron-electron interaction effects 0,0- Recently, two 
groups proposed that 2D arrays of superconducting is- 
lands immersed in a metal undergo a quantum supercon- 
ductor to metal transition (SMT) with decreasing island 
concentration |^. The resulting metallic phase has 
non-Fermi liquid properties including a pseudogap and 
anomalous magnetoresistance 0. 

Furthermore, theories of mesoscopic Superconductor- 
Normal metal (SN) structures have revealed mechanisms 
by which a finite DOS can appear within the energy gap 
of superconducting structures. These states, which give 
the DOS a hybrid superconductor-metal appearance, cor- 
respond to quasi-particles that become partially trapped 
in the N regions 'lO|, lUl Mesoscopic spatial fluctu- 
ations in the local conductivity can give rise to "quasi- 
localized" states within the N region. These states ap- 
pear within the gap, smearing the gap edge and creating 
a DOS down to Ep ^3- Within semi-classical models 
[nl I12I Us! , the DOS depends on whether the dynam- 
ics in the N region is chaotic or integrable. Quasiparti- 
cles on ballistic, integrable trajectories can become quasi- 
trapped in N regions and contribute a subgap DOS that 
grows linearly from Ep. 

We have conducted a series of experiments on ultrathin 
SN bilayers in an effort to observe the proposed SMT 

and its associated metallic phase. These bilayers 
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FIG. 1: Schematic cross section of the bilayer tunnel junction. 



can be driven toward the metallic phase by increasing the 
metal layer thickness, , at fixed superconductor thick- 
ness, ds- Within quasiclassical proximity effect theories, 
the bilayers ought to have a hard, BCS gap in the DOS 
and make a transition to a metal phase only as d^ — > 00 ■ 
Contrary to this expectation, our previous work showed 
that their superconducting transition temperatures, T^, 
decrease faster with d^ than quasi-classical predictions 
suggesting the approach to a SMT Q . 

Here we present electron tunneling measurements 
showing that the DOS of ultrathin SN bilayers develops 
a hybrid superconductor-metal appearance that becomes 
more metallic as Tc decreases. Specifically, the super- 
conducting gap of ultrathin Pb-Ag bilayers systemati- 
cally fills with states as Tc is decreased by increasing the 
metal thickness, dAg at fixed superconductor thickness, 
dpb. The subgap DOS is finite at Ep and rises linearly 
with energy with a slope that increases with dAg ■ Recent 
theories suggest that these subgap states are quasiparti- 
cles that become decoupled from the superconductor for 
times longer than the superconducting coherence time 
miil 12, 13]. The linear dependence is consistent with 
semiclassical theories of quasiparticles becoming trapped 
on integrable trajectories [lllll^ . 

For these studies, Pb/Ag bilayers with ultraclean 
Pb/Ag interfaces were fabricated and measured in situ 
using quench condensation techniques in the UHV en- 
vironment of a dilution refrigerator crvostat|l4j. The 
metals were thermally evaporated onto fire polished glass 
substrates held at 8 K. Au/Ge contact pads and oxidized 
Al counterelectrodes with a small amount of magnetic 
impurities to prevent them from superconducting were 
deposited prior to cryostat mounting. To form bilayers. 
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a thin < 6 nm, electrically discontinuous, Pb film was 
deposited first followed by, without breaking vacuum or 
warming, a series of Ag depositions (see Fig. 1). The 
latter drove the bilayer through an insulator to super- 
conductor transition 0, . This procedure yielded a 
series of bilayers with a single dpt and a range of c?Ag 
that were probed with the same, 1.25 mm^ area, tun- 
neling counterelectrode and barrier. The transport and 
tunneling measurements were performed using standard 
4 terminal, low frequency AC techniques. Data acquired 
on a series of bilayers with dpt =4.0 nm and 4.2 < dAg < 
19.3 nm are presented here. This data set is the most 
complete and systematic. It exhibits features that are 
similar to those of other series with 1.5 < dpb < 6.0 nm. 

The normahzed, superconducting bilayer DOS, Ns{E), 
is obtained from the normalized differential conductance, 
Gj, of the tunnel iunction[l7l|: 
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where Gs and Gn are the differential conductances, 
dl/dV, in the superconducting and normal states, re- 
spectively. / is the tunnel current and V is the volt- 
age across the junction. / is the Fermi function and 
E is the energy measured from Ep. At low tempera- 
tures^ < 0.1 Tc, w Ns{eV). In the Cooper limit 
El IllllSl, Ns is pred icted to as sume the BCS form, 
N^'^^{E,A) = Re{E/^{E^ ~ A2)) where A is the en- 
ergy gap. 

The evolution of the resistive transitions, R{T) and 
tunneling conductances, Gj, of a bilayer series driven to- 
ward the metallic state is shown in Fig. 2. The tran- 
sitions are sharp and Tc, defined as the temperature at 
which i?(T) is half its normal state value, drops exponen- 
tially with dAg (inset). The Gj, obtained at T = 60 mK 
■C Tc, qualitatively resemble the BCS form, exhibiting 
an energy gap structure consisting of symmetric peaks 
and a depression at low voltages (Fig. 2b). Within the 
gap region, however, Gj has an approximately linear volt- 
age dependence and a finite value at zero voltage (Fig. 
2c) rather than an exponentially small value. This sub- 
gap conductance grows with dAg and fills nearly 40% of 
the gap for the bilayer with Tc = 0.67 K. In addition, 
the conductance peaks are shorter and broader than the 
BCS prediction. 

Depositing Pb atop a bilayer to create a trilayer re- 
verses the above evolution as shown in Fig. 3 for a bilayer 
with dph = 1.4 nm and dAg = 7.1 nm. It's Gj resembled 
that of the lower Tc bilayers in Fig. 2. Adding an up- 
per Pb layer sharpened the peaks, increased the energy 
gap and reduced the slope and zero voltage bias value 
of Gj. A second Pb evaporation continued these trends. 
We hasten to note that the reduction in the subgap con- 
ductance induced by the upper Pb layer is a sign that 
the subgap conductance reflects an intrinsic feature of 
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FIG. 2: (a) Sheet resistance vs. temperature for Pb/Ag bi- 
layers with dpb = 4 nm and dAg ~ 4.2, 6.7, 9.1, 12.4, 15.6, 
19.3 nm. Inset: Semilog plot Tc vs. dAg- (b) Gj vs. V for 
the same bilayers. (c) Same data as in (b) on a finer voltage 
scale. The dashed lines are guides for the eye. 



the DOS. A subgap conductance stemming from leakage 
would be unaffected by an upper Pb layer. 

Except for the subgap conductance and the broadened 
conductance peaks, the data in Fig. 2 follow quasiclas- 
sical models of the proximity effect in the Cooper limit 
[lil nil I20I l2l|. In this limit, which applies to bilayers 
with dsjdjsr << ^, the superconducting coherence length. 
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FIG. 3: (a) Gj at T = 60 mK vs. voltage for a bilayer 
dpb/dAg = 1.4 nm/7.1 nm, and two trilayers, dpt/dAg/dpt, 
of 1.4 nm/7.1 nm/2.2 nm and 1.4 nm/7.1 nm/5.4 nm. (b) 
Same data as in (a) on finer scales. 



electrons near Ep pass back and forth between the S and 
N regions rapidly compared to A/Ti. The effective su- 
perconducting coupling constant, A oc ln(Tc) (x /n(A), 
is the volume average of the coupling constants in the N 
and S regions and thus, depends linearly on dAg in agree- 
ment with the inset of Fig. 2a j3| • The rapid motion 
renders the pairing amplitude uniform across the bilayer 
and thus, leads to a BCS form for the tunneling DOS [l9| 
that roughly agrees with the data. 

The extra breadth in the peaks can be attributed to 
the existence of a distribution of energy gaps in the bi- 
layers. A distribution can naturally arise due to spa- 
tial variations in the film thickness ratio x — dAg/dpb 
since A oc Tc depends exponentially on x. Presuming 
A = Aoexp{—kx), a (random) normal distribution of x 
yields a log-normal distribution for A and a broadened 
form of Ns{E): 
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where A is the most probable energy gap and a is the 
width of the x distribution. The parameters Aq =0.88 



meV and k = 0.46 were obtained by estimating A for 
each bilayer as the voltage at which Gj =1 and fitting to 
the exponential form. Calculated Gj (see Fig. 4a) with 
broadened peaks that resemble the data require a <25% 
of X. 

The states that systematically fill the superconducting 
gap, however, do not readily conform to a gap distribu- 
tion model. The log-normal gap distribution drops too 
rapidly below A to reproduce the linear DOS and finite 
zero bias conductance(see Fig. 4a). Creating the ob- 
served DOS at Ep would require that a assume values 
at least 10 times larger than the values used to fit the 
peaks. Thus, we are led to the conclusion that the sub- 
gap DOS corresponds to quasiparticles that fall outside 
the Cooper limit picture and thus must relate to a new 
physical mechanism. A similar problem was encountered 
by Gupta and coworkers |22j who were unable to fit sub- 
gap structure in their STM data on Au/Nb films. 

According to semi-classical theories, subgap states cor- 
respond to quasiparticles that propagate in an N region 
over distances more than ^ between successive Andreev 
scattering events with an NS interfaceHl [H IH 
This situation normally occurs in SN structures with N 
region dimensions exceeding ^, |3 li^l EH l2^ or 
structures with integrable dynamics Il2l |. in which 
quasiparticles can execute nearly closed trajectories. In- 
terestingly, the predicted DOS grows approximately lin- 
early with energy for a variety of structures |TlL Il2. HH 
In particular, Melsen and coworkers |llj proposed 
that for a rectangular integrable N billiard attached 
to a bulk S region, the distribution of path lengths, 
P{L), between successive Andreev scattering events fol- 
lows P{L) (x a,s L ^ oo. The normalized slope 
of the resuhing linear DOS is 2/{-KETh) im, where Eph 
is the Thouless energy. Furthermore, the more general 
case of a mixed phase space that consists of both regu- 
lar (integrable) and chaotic regions also yields a linear 
DOS and, in addition, a constant background DOS aris- 
ing from phase space regions that arc completely discon- 
nected from the superconductor (l.H] . 

Guided by the above considerations and low tempera- 
ture STM measurements showing that the bilayers con- 
sist of crystalline Pb and Ag grains |2^ [s^, we now as- 
sume that the bilayers present a mixed phase space of reg- 
ular and chaotic regionsj36j. Quasiparticles with chaotic 
trajectories contribute a BCS like portion, iV^, to the 
DOS while those with regular trajectories contribute a 
linear subgap DOS. Thus, we parameterize Ns{E): 



NsiE) 



N^{E,A) ifE>Ec 
aE + l3 if E <Ec 



where a and f3 are the slope and the intercept of the 
linear dependence, respectively, and Ec is the energy at 
which the two DOS intersect. Fig. 4a(right) shows Gj 
calculated from this DOS Fits to the data opti- 

mized to capture the subgap slopes and the peak heights 
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are shown in Fig. 4b and Fig. 4c. A decreases expo- 
nentially with dAg (inset Fig. 4b) consistent with expec- 
tations for the Cooper limit. Since a is the slope of the 
DOS averaged over the regular phase space regions, a is 
proportional to the average E^]^. For these bilayers, we 



expect 



d\g/{hD), where D is the intergrain dif- 



fusivity and thus a oc d\g. Fig. 4c compares the exper- 
imental a values with those estimated by 2d\g / {ttTlD) ^ 
using _D = 5 X l{)~^m? as obtained from transport 
measurements. The two roughly agree at small dAg and 
deviate at higher dAg where the experimental a values 
grow more rapidly than expected. The deviation sug- 
gests that a also depends on Tj.. This case has not been 
considered by the current theories 0, . The trilayer 
data support this conjecture, exhibiting a reduction in a 
when Tc is increased at fixed dAg- Regardless, this devia- 
tion implies that the regular phase space regions occupy 
an increasing fraction of the total phase space volume as 
dAg increases. The concomitant increase of (3 implies a 
growing fraction of disconnected regions, as well. Conse- 
quently, as Tc decreases, a growing fraction of quasipar- 
ticles decouple from the S layer for times much greater 
than the superconducting coherence time. 

Qualitatively, the growth in the subgap DOS gives 
the DOS a hybrid-metal-superconductor appearance that 
may be a sign of an approaching SMT. Ultrathin films 
near the disorder tuned superconductor-insulator transi- 
tion exhibit a similar filling of the gap 113 and low 
energy states have been invoked as a source of dissipa- 
tion that drives the Quantum SMT observed in nanowires 
[33I I34L 35.J . Perhaps these states arise due to similar 
quasiparticle trapping effects. 

Tunneling experiments on ultrathin Pb-Ag bilayers at 
low reduced temperatures (T/Tc < 0.1) have revealed 
an unexpected linear DOS within their superconducting 
energy gap. The fraction of the gap filled with states 
grows with increasing dAg- We have identified the sub- 
gap states as quasiparticles that are weakly coupled to 
the superconductor layer. Within a semiclassical picture, 
these quasiparticles become trapped in regular regions of 
phase space. 

The work was supported bv Inst-dmr020360'8l We ac- 
knowledge helpful conversations with Dmitri Feldman. 



FIG. 4: (a) Comparison of Gj calculated using a log-normal 
gap distribution (left) and a log-normal distribution combined 
with a linear DOS within the gap (right). The corresponding 
curves have the same fitting parameters, A and a. (b) Data 

from Fig. 2a with fits. The fitting parameter o"/(^^) de- 
creases from 0.23 to 0.11 as cLaq goes from 4.2 nm to 19.3 nm, 
while 13 increases from 0.033 to 0.065. Inset; Semilog plot of 
A, vs. dAg/dpb- (c) Same as (b) on finer scales. Inset: a vs. 
d'Ag. The solid line gives the prediction based on transport 
measurements. 
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